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In a pioneer experiment, Bohlcin et al. realized the controlled sliding 
of two-dimensional colloidal crystals over laser-generated periodic 
or quasi-periodic potentials. Here we present realistic simulations 
and arguments which besides reproducing the main experimentally 
observed features, give a first theoretical demonstration of the po- 
tential impact of colloid sliding in nanotribology. The free motion 
of solitons and antisolitons in the sliding of hard incommensurate 
crystals is contrasted with the soliton-antisoliton pair nucleation at 
the large static friction threshold F s when the two lattices are com- 
mensurate and pinned. The frictional work directly extracted from 
particles' velocities can be analysed as a function of classic tribo- 
logical parameters, including speed, spacing and amplitude of the 
periodic potential (representing respectively the mismatch of the 
sliding interface, and the corrugation, or "load"). These and other 
features suggestive of further experiments and insights promote col- 
loid sliding to a novel friction study instrument. 
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The intimate understanding of sliding friction, a cen- 
tral player in the physics and technology of an enor- 
mous variety of systems, from nanotribology to mesoscale and 
macroscale sliding [UH], is historically hampered by a num- 
ber of difficulties. One of them is the practical inaccessibility 
of the buried interface between the moving bodies - with few 
exceptions, we can only hypothesize about its nature and be- 
havior during sliding. Another is the general impossibility to 
fully control the detailed nature, morphology, and geometric 
parameters of the sliders; thus for example, even perfectly 
periodic, defect-free contacting surfaces have essentially only 
been accessible theoretically. If we knew and, on top of that, 
if we could control the properties and the relative asperity 
parameters of the sliders, our physical understanding could 
greatly increase, also disclosing possibilities to tune friction 
in nano and mesoscopic systems and devices. As Bohlein et 
al. [3] showed, two dimensional (2D) colloid crystalline mono- 
layers can be forced by the flow of their embedding fluid to 
slide against a laser-generated static potential mimicking the 
interface "corrugation" potential in ordinary sliding friction. 
The external pushing force, the interparticle interactions, and 
especially the corrugation potential are all under control, the 
latter ranging from weak to strong, and from periodic, to 
quasi-periodic [4] [5], in principle to more complex types too. 
Contrary to established techniques in meso and nanosize slid- 
ing friction (Atomic Force Microscope, Surface Force Appa- 
ratus, Quartz Crystal Microbalance) [BJ, which address the 
tribological response in terms of averaged physical quantities 
(overall static and kinetic friction, mean velocities, slip lengths 
and slip times, etc.), in colloid sliding every individual particle 
can in principle be followed in real time, stealing a privilege 
hitherto restricted to the ideal world of molecular-dynamics 
(MD) simulations 0[U E|. 

Materializing concepts long-anticipated theoretically |10l 
111] , the colloid sliding data showed how the sliding of a flat 
crystalline lattice on a perfectly periodic substrate takes place 
through the motion of soliton or antisoliton superstructures 



(also known in one dimension - ID - as kinks or antikinks) - 
positive or negative density modulations that reflect the mis- 
fit dislocations of the two lattices that are incommensurate 
in their mutual registry. While forming regular static Moire 
superstructure patterns when at rest, solitons constitute the 
actual mobile entities during depinning and sliding, and are 
essential for "superlubricity" [12] - i.e., zero static friction - 
of hard incommensurate sliders. When solitons are absent 
at rest owing to commensurability of the two sliders (or are 
present but pinned in soft incommensurability), the colloids 
and the periodic potential are initially stuck together. Only 
after the static friction force F s is overcome, solitons appear 
(or depin if they already exist but are pinned) unlocking the 
colloids away from the corrugation potential, so that sliding 
can take place. 

Our aim here is to understand and demonstrate, based 
on molecular dynamics (MD) sliding simulations, how the 
great colloid visibility and controllability can be put to di- 
rect use in a tribological context. The full phase diagram 
versus colloid density and sliding force is explored first of 
all, highlighting a large asymmetry between solitons and 
antisolitons, and a strong evolution from commensurate to 
incommensurate caused by sliding. We then extract and 
predict the frictional work as a function of mean veloc- 
ity and corrugation amplitude, loosely mimicking "load" - 
the same variables of classic macroscopic friction laws. We 
also discuss new local phenomena underlying depinning, in- 
cluding the edge-originated spawning of incommensurate an- 
tisolitons and the bulk-originated nucleation/separation of 
soliton-antisoliton pairs in the commensurate well as 

global analogies to driven Josephson junctions, charge-density 
waves, and the sliding of adsorbate islands on crystal surfaces 
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Modeling and simulations 

The driven colloids are modelled as charged point particles un- 
dergoing overdamped 2D planar dynamics under an external 
force F, parallel to the plane, applied to each colloid. While 
the fluid is not described explicitly, F is to be interpreted as 
r/Vd where rj and Vd are the effective fluid viscosity and veloc- 
ity. Particles repel each other with a screened Coulomb in- 
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Fig. 1. (a) A sketch of the model for the colloid particles interacting with a 
periodic potential W . (b) The static initial configuration at p = 0.95. Three 
families of antisoliton lines (darker areas) cross at 120 degrees, (c) Velocity-force 
characteristics for various colloid densities, with a lattice-potential corrugation com- 
mensurability ratio p = ai as /a C oll = 1-0 (CO), 1.05 (SI), and 0.95 (Al). The 
CO case always displays static friction. For weak corrugation (Uq = 0.1), F s = 
in both Al and SI incommensurate cases. At larger corrugation (t/fj = 0.5) a major 
asymmetry appears between the Al and SI configurations: only the Al case exhibits 
a finite depinning threshold with static friction. (d,e) Snapshots of the central region 
of the initially commensurate colloid during motion, illustrating sliding-generated soli- 
tons, whose density increases as F is increased. In all snapshots, colloids located at 
repulsive spots of the corrugation potential [defined by W^r) > —Uq/2, e.g. the 
colloid pointed at by the red filled arrow in panel (a)] are drawn as dark red spots, 
while colloids nearer to potential minima [H^(r) < —Uq/2, e.g. the colloid pointed 
at by the blue empty arrow in panel (a)] are light blue. 
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Fig. 2. 2D Aubry transition for antisolitons at p = 0.95, in an infinite-size 
colloid system. Main panel: colloid mobility as a function of the applied force, for 
increasing corrugation amplitude Uq. Note the appearance of pinning with static 
friction just above Uq= 0.2. Inset: static friction (depinning) force F s , normalized 
to the single-colloid force barrier F s i , as a function of Uq, with an arrow indicating 
the critical Aubry corrugation. 

set this balance so that the average colloid lattice spacing a co ii 
(before submittal to the corrugation potential W) is unityQ 
We then realize a variety of mismatched ratios p — ai a a/a C oii 
by changing the corrugation period ai as . In the following, 
we focus on three representative cases, namely: underdense, 
p = 0.95 (antisoliton-incommensurate - Al; the starting state 
at rest is shown in Fig. Hp); ideally dense, p = 1.0 (nearly 
commensurate - CO, which becomes exactly commensurate 
after turningon W); overdense, p = 1.05 (soliton incommen- 
surate - SI) Fj. In order to simulate experiment, and also to 
prevent solitons from leaving the finite-size sample, our exter- 
nal force F is ramped in time in small well-spaced steps, so 
that its overall value slowly alternates in sign, forth and back 
with a long time period. Full simulation details are given in 
the supporting information. 



terparticle repulsion, V(ry) = Q/r,j exp(— ry/Au) with Ad 
substantially smaller than the mean distance between par- 
ticles. Colloids are immersed in a Gaussian-shaped over- 
all confining potential G(jrj) = — A c exp(— r' 2 /a 2 ) - the 
large radius a representing the laser spot size - and in a 
triangular-lattice periodic potential W(r) = — (2C/o/9)[3/2 + 
2 cos {2irr x I ai as ) cos(27rr M / (v^Oias ) ) +cos(47rr y / (villas ) )] rep- 
resenting the interface "corrugation", see Fig. Finally, in 
addition to the external force, a Stokes viscous force —rjWi 
acts on each particle i = 1, N, and accounts for the dis- 
sipation of the colloids kinetic energy into the thermal bath. 
We typically simulate N ~ 30, 000 - a particle number much 
smaller than in experiment, but sufficient to extract reliable 
physical resultsQ In the absence of corrugation (Uo = 0), col- 
loids form a 2D crystalline island at rest. The 2D density of 
the triangular 2D lattice is fixed by N and by the balance of 
the confining energy G and the 2-body repulsion energy. We 




Fig. 3. Depinned, moving particles (darker) of a 2D colloid in a periodic potential 
under the action of a rightward force F . (a) Rightward propagating solitons of over 
dense colloids (p = 1.05, SI); (b) Leftward propagating antisolitons of underdense 
colloid (p = 0.95, Al). 



To reduce simulation sizes and times, our particles form an island near the center of the Gaussian 
potential, which is the region experimentally visualized. Particles outside this region, whose role 
is less relevant, are omitted. Moreover, thermal effects (although straightforward to introduce in 
simulation), have not been accessed in colloid experiments. After verifying that the main features 
are not washed out at 300 K (see supporting information), here we will, for the sake of clarity, only 
present results obtained with a dissipative T — Langevin dynamics. 

2 The spacing of the fully relaxed colloid configuration varies smoothly from a ~ 0.984 at the 
sample center to a 1 .05 at the side, with an average density equal to that of a triangular crystal 
of spacing a co n — 1. 

Similar models were studied in the past with a view to understand two dimensional Frenkel- 
Kontorova models and adsorbate monolayers physics "20 21 22 23" 
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Results 

Figure QJ displays the mean speed (v cnx ) of the central por- 
tion of the colloid system as a function of the driving force 
F. Fully reproducing experiment [3], the simulated force- 
velocity characteristics of Fig. [TJ; show a large static fric- 
tion force threshold in the p ~ 1 CO case, where the colloid 
and corrugation lattices are pinned together. Static friction 
is lost in case of incommensurability and moderate corru- 
gation, where preformed mobile solitons or antisolitons are 
present. The snapshots of Fig. [3] illustrate the patterns of 
solitons/antisolitons sliding in opposite directions under the 
same driving force F > 0. For a weak external force and a 
~ 5% lattice mismatch, the static friction drops essentially to 
zero, and a nearly free viscous sliding is realized, reflecting 
a situation of "superlubricity" [241 1251 126| . However, under 
the same conditions, not all incommensurate geometries are 
superlubric. Whereas for weak corrugation the overall colloid 
mobility (v cul ) /F is remarkably constant for both incommen- 
surate densities, we find in fact that by increasing the cor- 
rugation amplitude Uq the mobility of the AI configuration 
drops to zero at small force, and pinning with static friction 
reemerges despite incommensurability. By contrast, SI con- 
figurations remain superlubric up to much larger Uo- 

The Aubry transition. Borrowing results of the ID Frenkel 
Kontorova (FK) model |10] . the single-soliton width d ~ 
g 1/2 a la _ s where g = £af as fc/(7o [here k = V"(a coU ) and £ is 
a constant of order unity] is large for a hard layer on a weak 
corrugation, and small for a soft layer on a strong corruga- 
tion. Between these two extremes, the ID incommensurate 
FK model crosses the so-called Aubry transition [24] where 
superlubricity is lost, and pinning sets in with static friction 
despite incommensurability. Even in the present 2D case it is 
qualitatively expected that all incommensurate colloids, both 



underdense (p ^ I) and overdense (p ^ 1) will undergo an 
Aubry-like superlubric-to-pinned transition for increasing cor- 
rugation. 

This expectation is confirmed in our 2D model colloid 
system. Figure [2] (obtained by independent simulations of 
the infinite-size system with periodic boundary conditions) 
shows the Aubry-like pinning transition crossed by an AI 
(p — 0.95) underdense colloid at a critical corrugation, here 
U(j rlt ~ 0.2 — 0.3. The threshold Aubry corrugation depends 
upon p, and is much larger for overdense SI than for under- 
dense AI colloids. 



Soliton-antisoliton asymmetry.This strong asymmetry of 
static friction - and of all other properties - between overdense 
(.P lp 1) an d underdense (p < 1) colloids can be rationalized, 
in the limit of strong corrugation g <C 1, in terms of the large 
physical difference between solitons, defects formed by lines of 
lattice interstitials, and antisolitons, lines of vacancies. This 
asymmetry remains even for weak corrugation (j> 1), when 
solitons/antisolitons involve relative displacements far smaller 
than those of proper interstitials or vacancies. A small varia- 
tion 8 in the inter-colloid separation a is sufficient to produce 
a large relative variation of the effective spring constant, i.e. 
the interaction curvature 
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For a realistic Ad = 0.03 a co n, this highly nonlinear and asym- 
metric relation, implies a huge 460% increase whenever two 
colloids are approached by 5% of their average separation, 
but only a 82% reduction for a 5% increased separation. This 
asymmetry is held responsible for the much weaker propen- 
sity of solitons to become pinned and to localize compared to 
antisolitons. 



The sliding state. Under sliding, the shapes and geometries 
of solitons/antisolitons and their motion are of most immedi- 
ate interest, as they are directly comparable with experiment. 
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Fig. 4. Three successive snapshots of the initial depinning instants of the com- 
mensurate (p = 1) configuration for the F C± F s i simulation, see Fig. The 
horizontally extended window visualizes the the nucleation and separation of a soliton- 
antisoliton pair, left of the central observation region (square). Pair nucleation con- 
stitutes the depinning mechanism of all commensurate sliders. 



Fig. 5. Effective "phase diagram" of the finite colloid crystal as a function of 
the lattice spacing mismatch p, for Uo = 0.1. The two potential energy curves 
(colloid-colloid repulsion plus W interaction) characterize the static phases: (red 
squares) relaxation started from the slightly inhomogeneous configuration produced 
by a previous relaxation for Uq = 0; (blue circles) relaxation started with a fully 
matched lattice of spacing Gti as . Two AI (underdense) and SI (overdense) phases 
surround a commensurate phase (CO). The depinning mechanisms and the soliton 
structures sustaining sliding are illustrated at top of figure for the different regions. 
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Figure [3] shows the large-scale checkerboard structure of soli- 
tons/antisolitons of the sliding colloid lattice. They move 
with a speed v much larger than the average lattice speed 
(t>cm), because v/{vcm) ~ p/{p ~ 1) by particle conservation. 
The moving structure is a distortion of the original triangular 
soliton/antisoliton pattern (Fig. QJ)) induced by the circular 
shape of the confining potential, and by the directional slid- 
ing. With increasing F, the soliton arrangements elongate 
into a stripe-like pattern perpendicular to the driving direc- 
tion. Comparison with experimental pictures is quite realistic, 
especially when focusing (as done in experiment) on the cen- 
tral sample region, far from boundaries. 

In the AI superlubric colloid p < 1, preformed antisoli- 
tons fly (leftward) across the colloid lattice antiparallel to the 
(rightward) force. They are eventually absorbed at the left 
edge boundary, while new ones spawn at the right edge bound- 
ary to replace them, sustaining a steady-state mobility. In the 
SI superlubric colloid p > 1 conversely, preformed solitons fly 
rightward, parallel to the force. Solitons, unlike antisolitons, 
are not automatically spawned at the boundary, owing to the 
decreasing density. Instead, an antisoliton/soliton pairs must 
nucleate first, near the boundary, and this is possible only if 
the force overcomes the nucleation barrier. Below this thresh- 
old, we observe that a steady DC external force eventually 
sweeps out all the preformed solitons transforming the colloid 
to an artificially pinned, immobile CO state. 

Finally, the pinned CO colloid p > 1 only moves after 
static friction is overcome. As illustrated in Fig. [4l motion 
starts off here by nucleation of soliton-antisoliton pairs inside 
the bulk - here close to the left edge because the central region 
tends to be slightly overdense. The antisolitons flow leftwards 
and are absorbed by the left edge, becoming undetectable to 
the optically monitored central part of the colloid, where only 
solitons transit, as seen in experiment. This type of commen- 
surate nucleation has been described in considerable detail 
in literature, including finite-temperature effects [111 127] . We 
note here that in the pinned CO colloid the soliton or antisoli- 
ton density, initially zero, actually increases with increasing 
sliding velocity (see, e.g., Fig. [ljl,e), as opposed to frankly 
incommensurate cases, where it is nearly constant. 



In the running state, the colloid average density increases or 
decreases from 1 to a value closer to the nominal p of the col- 
loid at Uq — 0. This explains why in a quasi-commensurate 
configuration with p > 1 such as that shown by Bohlein el 
al. [3], solitons (and not, e.g., soliton-antisoliton pairs) sweep 
the colloid upon depinning, as also seen in Fig. [TJi,e. 

It is curious to note here the different fate of solitons in 
the slightly overdense CO and in the SI phases. In the CO 
phase they do not exist at rest, but they appear after depin- 
ning and under sliding. In the SI phase they exist at rest, 
but they could be swept out under DC sliding, when a weak 
external force can turn the SI colloid into effectively CO. We 
never saw this sweepout phenomenon on the AI side. 

Frictional analysis 

We turn now to frictional work, a quantity of crucial impor- 
tance for the tribological significance of colloid sliding. We 
can write the overall power balance as the scalar product of 
the instantaneous velocity v, of each colloid i by the net force 
acting on it, r\ (vj — Vi), including both the bare external force 
F = r\\d and the viscous drag —rjVi. This product vanishes 
instantaneously at any time when either colloids are stuck 
(vj = 0) or else when the corrugation potential is absent, so 
that Vi = Vd. After averaging over a very long trajectory, the 
balance reads 

ftot = ^{(Vd - Vi) ■ Vi) 

i 

= (JVF ■ <v cm ) - r,(|v cm | 2 )) - r^du,! 2 ) 

i 

= Pfrict-Pkin, [2] 

where u; = Vj — v cm . Under steady-state sliding conditions 
where Ptot = 0, the effective friction power Pf r i ct is exactly 
balanced by an internal kinetic energy excess rate. Per col- 
loid particle, Pf r i ct is 

Pwct = 2i F ■ (v cm ) - r?|(v cm )| 2 , [3] 



Phase-diagram evolution with sliding. Much can be learned 
about the habit of sliding colloids from their behavior and 
their structural phase diagram, first at rest and then under 
sliding. With p ~ 1, close to commensurate but not ex- 
actly commensurate, the colloid monolayer can realize in the 
periodic potential two alternative static arrangements which 
are local minima of the overall free energy: a fully lattice- 
matched CO state, or a weakly incommensurate state charac- 
terized by a sparse soliton (AI or SI) superstructure, with a 
density fluctuating around the local value prescribed by the 
G — V balance. Comparing the potential energy of these two 
states as a function of p, the static phase diagram contains, 
as sketched in Fig. O a fully commensurate extended CO re- 
gion separated from the AI and SI regions by commensurate- 
incommensurate transitions, well known in adsorbed surface 
layers [28] [29j [30l [31] . The CO region is wider on the SI side 
(p > 1) than the AI side (p < 1), another manifestation of 
the SI-AI asymmetry discussed above. The CO range natu- 
rally widens or shrinks when the corrugation amplitude Uo is 
increased or decreased. 

Under an external force F, sliding effectively tilts the bal- 
ance between the two static phases (loosely speaking, for of 
course under sliding the physical significance of a "phase" is 
not the same as at rest) sliding populates the former CO phase 
with solitons/antisolitons, turning it effectively into SI or AI. 
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Fig. 6. Sliding friction Pf r i c t per particle as a function of speed {v C m} for an 
underdense AI colloid p = 0.95 for increasing corrugation amplitude Uq. Inset: Uq 
dependence of friction for speed (v cm ) = 0.3, showing the quadratic rise for weak 
corrugation behavior, followed by a roughly linear growth. 
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where small center of mass fluctuations are neglected, by as- 
suming (v cm 2 ) ~ |(v cm )| 2 . 

Figure [6] shows Pf r i ct (briefly referred to as "friction" in 
the following) for the special case of the AI underdense phase 
as extracted as a function of (v cm ) through a "bulk" simula- 
tion (with periodic boundary conditions as in Fig. [2}. The 
main features found are (i) a linear rise at low CM speed; (ii) 
a decline at large speed; (iii) a maximum at some intermedi- 
ate corrugation-dependent speed. We moreover observe that 
(iv) the dissipated power increases (not unexpectedly) with 
corrugation; and (v) the corresponding frictional maximum 
simultaneously shifts to larger speed. 

The qualitative interpretation of these results is relatively 
straightforward, and yet revealing, (i, iv) At low sliding ve- 
locities the motion of solitons/antisolitons involves the vis- 
cous motion of individual particles with a velocity distribu- 
tion whose spread toward higher values rises proportionally 
to the sliding speed and inversely proportional to their spatial 
width. As shown, e.g. within the ID FK model [10], but also 
in the present simulations, the width d of solitons/antisolitons 
increases roughly as d ~ a^^g with the dimensionless inter- 
particle interaction strength g oc a 2 as V"(a co ii)/C/o measured 
relative to the periodic corrugation amplitude. The decrease 
in width with increasing corrugation Uo requires an increas- 
ing instantaneous speed of individual particles in the soli- 
ton/antisoliton, yielding an increasing viscous friction, and 
a decreasing overall mobility as observed, (ii, iii, v) At high 
sliding velocities, the colloid relaxation time exceeds the soli- 
ton/antisoliton transit time across the Peierls-Nabarro barrier 
[10| so that their spatial structure is gradually washed out by 
the sliding motion. The critical speed where the smoothen- 
ing behavior takes over, roughly corresponding to maximal 
friction, increases as corrugation increases, corresponding to 
narrower solitons/antisolitons that are harder to wash out. 
The increase of friction with corrugation strength Uo, plot- 
ted in the inset for a chosen speed, is found to be quadratic at 
weak corrugation, gradually turning to linear for larger values. 
Linear response theory naturally accounts for the quadratic 
increase, a behavior first discussed by Cieplak el al. [32] and 
observed in quartz crystal microbalance experiments [33] . 

Demonstrated for a specific AI case with antisolitons, the 
above results appear in fact of general validity for infinitely ex- 
tended sliders of controlled colloid density, and apply equally 
well although with great quantitative asymmetry to SI with 
solitons, once their larger widths and greater mobilities and 
weaker Peierls-Nabarro barriers are taken into account. 



Summary and Conclusions 

In this study we have presented initial simulation results and 
theory that strongly vouch in favor of sliding of colloid layers 
on laser-originated corrugations as a promising tool for future 
tribological advances. The motion of solitons and antisolitons 
known from experiment is reproduced and understood, un- 
raveling the subtle depinning mechanisms at play. The pres- 
ence of Aubry transitions is pointed out for future verification, 
along with a strong asymmetry between underdense and over- 
dense incommensurate layers. Of direct tribological interest, 
we anticipate the behavior of friction with corrugation (mim- 
icking "load") and with sliding velocity, with results which, 
while of course generally very different from the classic laws 
of macroscopic friction, are highly relevant to friction at nano 
and mesoscopic scales. Our approach moreover indicates a 
strong complementarity between theory plus simulation, and 
experiment, an aspect which we intend to pursue further. 

There are many lines of future research that this study 
implicitly suggests. One line will be to pursue the analogy of 
the sliding over a periodic potential with other systems such 
as driven Josephson junctions [34] , and sliding charge-density 
waves [2]. Time-dependent nonlinear phenomena such as 
the Shapiro steps [141 134] should become accessible to colloid 
sliding too. A second line is to include non-periodic complica- 
tions to the corrugation potential, including the quasicrystal 
geometry such as that recently realized [35] and beyond that, 
random, or pseudo-random corrugations to be realized in the 
future. A third line involves the investigation of the lubri- 
cant speed quantization phenomena, characterized so far only 
theoretically [36l [371 [381139] . 

A further very important development will be to address 
colloidal friction in larger, mesoscopic or macroscopic size sys- 
tems, whose phenomenology is accessible so far only by a few, 
very ingenious, but very limited, methods [6] 1161 1171 1401 1411 
142] . A major scope in that case will be to realize and study 
stick-slip friction and aging phenomena, at the heart of real- 
istic physical and technological tribology. 
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